Dark Excitons in Carbon Nanotubes 
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We investigate the lowest many-body excited states in carbon nanotubes by means of ab initio 
calculations. On the basis of these calculations and an additional theoretical analysis of the excitons, 
we demonstrate that the splitting between the dark and bright spin-singlet states is due to the direct 
Coulomb electron-hole interaction, that the dark exciton always has gerade symmetry, and that it is 
always lower in energy relative to the bright exciton. Applying these studies to 0.4-0.8-nm-diameter 
tubes, we find a gerade-ungerade splitting of several tens of an meV, and a singlet-triplet splitting 
of comparable value. Our symmetry analysis of the origin of the singlet splitting indicates possible 
startegies to increase the emission efficiency of nanotubes. 
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In recent years, the luminescence properties of carbon 
nanotubes have sustained intense theoretical 0,0,13 and 
experimental research 0, IS IS 0] Of particular 
interest is the study of the emission efficiency and the 
radiative lifetimes, in view of applications in the field of 
opto-electronics and photonics 0, 0, El 03 ■ 

Recent studies on excited-state dynamics Hij . 
combined with the experimentally observed low quantum 
efficiency, indicate that exciton relaxation is dominated 
by non-radiative processes, which dramatically quench 
emission spectra. The presence of optically forbidden 
states, i.e., dark states, below the optical gap, to which 
the system can decay non-radiatively, is believed to be 
responsible for the luminescence quenching 0, Hr| . 
Their presence is estimated to increase the effective ra- 
diative lifetime of about five times ^S]- However, evi- 
dence of dark states has been presented only recently in 
experiments involving large magnetic fields |17j , and very 
little is known about this subject. 

Thus far, no theoretical mechanism has been proposed 
explaining the reason why there are dark states lying 
below the first bright one nor if this is a general rule for 
all kinds of tubes. 

In this work, we carry out ab initio calculations of ex- 
citonic states for several tubes of sizes and dimensions of 
interest in the experimental physics community. On the 
basis of our calculations, we show that only two pairs of 
bands are involved in the description of the lowest exci- 
tons. This fact, as we shall see later, implies that there 
is always a dark singlet state below the lowest bright 
one. Finally, our symmetry analysis allows us to pro- 
pose alternative ways to increase the quantum efficiency 
of nanotubes through symmetry- breaking mechanisms. 

The excited states are calculated using a many-body 
theoretical approach which involves both self-energy cor- 
rections in the treatment of single-particle energies and 
effects arising from the electron-hole interaction. The 
former is treated within the GW approximation [lsl llflj 
while the latter is obtained from the diagonalization of 
the Bethe-Salpeter Hamiltonian [5( 1 EH U2. . 





FIG. f : The two principal symmetries of carbon nanotubes: 
the roto-translation, RT, with axis z, and the 7r-rotation 
about the axis U, Ru- Left panel: RT is shown by the dark 
atoms. The axis U pierces the center of an hexagon and is per- 
pendicular to the screw axis z. Right panel: The two inequiv- 
alent atoms A and B of each hexagon are shown in greater 
detail. When Ru is applied to the atomic positions, A goes 
into B and viceversa. 



In our calculations, we take advantage of all the crystal 
symmetries of the tube (see Fig.^l, which give us a better 
understanding of the band structure as well as of the 
excited states 0j US ■ A typical band structure of the 
tubes, shown in the upper panel of Fig. [21 in the extended 
zone scheme, has a minimum gap at the points labelled 
K and K'. 

In all the tubes studied, for an energy window from 
to 5 eV, which is the range of interest for optical ab- 
sorption experiments, the only transitions that come into 
play for the lowest singlet states are those between the 
valence band maximum and the conduction band mini- 
mum for fc-points near K and K' . In fact, for the singlet 
states lit and lg, where "1" refers to the lowest bound 
excitonic state and u (g) refers to the parity with respect 
to 7r- rotations about the axis U (Fig. Q), most of the 
weight (sa 99%) is concentrated near the gap in K and 
K' . It can be clearly seen in Fig.0 where the probability 
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FIG. 2: A typical band structure (upper panel) with the 
wavefunctions for the two lowest excitons (two lower panels) 
in the extended zone scheme. Both wavefunctions are local- 
ized around K and K' and manifest even and odd parity with 
respect to the Ru symmetry. The diagrams shown are for the 
(6,4) tube. 
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FIG. 3: The lowest singlet states in the absorption spec- 
trum of the (6,4) tube. The shaded line is the single-particle 
absorption, ignoring excitonic effects. The solid and dotted 
vertical lines indicate the positions of the two lowest excitons 
(Ig and 1m), where the lower of the two is dipole forbidden. 



amplitude of the transition is plotted as a function of k 
in the extended zone scheme. In addition, we find that 
the dark states have the same excitonic composition, as 
far as the probability amplitude is concerned, since the 
bright states differ from them only in the phase informa- 
tion (lowest panel in Fig.|2). This behaviour suggests the 
possibility of describing the lowest excitonic states with 
a graphite-derived, four-band model. 

Let us call the excitonic wavefunction concentrated 
near K, Fk, and the corresponding wavefunction at K\ 
Gfe. Let the indices i and j denote the space of possible 
optical transitions vk — > ck and H^, the Bethc-Salpcter 
(BS) Hamiltonian in this space of transitions. Fk and Gk 



can be written alternatively in the basis of the transitions 
as Fi and Gj. 

Hij can be decomposed in three parts, 

Hij = Dij + 2Xij — Wij = Dij + Fj . 

Here is the diagonal term, X+j is the exchange term, 
and Wij is the direct Coulomb term. Fj = 2Xij — Wij is 
the total electron-hole interaction term. 

We define the parameters E = E gap — Eluding and A, 

as 

E = Y J FtH ij Fj = Y,G*H ij Gj (1) 



and 



& = Y J F*H ij G j = Y,F*hGj. (2) 



We can obtain excitonic states diagonalizing the matrix 

It has two eigenvectors, v±, given by 

v± = 71 ( ±A*/|A| 

with eigenvalues, A± = E ± |A|. 

We calculate A± using our ab initio BS hamiltonian 
and model wavefunctions, which take into account only 
contributions from transitions at k-points very close to 
K and K'. The splitting | A| is found with an error of no 
more than 1% with respect to the ab initio calculation of 
the singlet splitting Ef u — Ef . We can therefore safely 
conclude that no bands other than the four graphite- 
derived ones come into play for the lowest singlet exciton 
states and the model assumption is correct. 

We have yet to demonstrate which of these two states, 
v±, is bright. In order to simplify the following discus- 
sion, note that the direct term W is large compared to 
the exchange term X, which will therefore be neglected 
in what follows with an error less than 1 %. 

We define the real-space excitonic wavefunction 

* K {v,r') = Y,F k ct>* ck {r)ct> vk {r'), 

k 

where </>* fc (r) and (/» 1) fe(r / ) are the single-particle conduc- 
tion and valence wavefunctions. 

Without loss of generality, we examine only the term 
with k — K, therefore 

vMr,r')~^ fe (r)<Mr')- 

We exclude the modulation factor Fk for simplicity, since 
the argument below would still be valid. Since the sym- 
metry Rjj depicted in Fig. ^ transforms K in K' , and 
consequently, Fk in Gk-, then Eq. [21 can be written as 



A = 



J d\ J d 3 r'y* K (r,r')W(r,r')y K {~Rur,-R v r'). 
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If we assume that 4> c k and <$> v k are 7r-like, i.e., are 
symmetric (7r) and antisymmetric (tt*) LCAO combina- 
tions of basis functions centered on the two inequiva- 
lent atomic sites A and B (which we denote as Xa{ y ) 
and XB{r)), then cA -(r) = ^[xaW - e'^eW] and 
^(r) = ^b(r) + e l5 XB (r)]. 

From the relations xa(RuI') = Xs( r ) and Xb(RuI") = 
Xa( t )> we have 

^-(Ru^Rur') = -*^(r,r'). 
Therefore A can be written as 

A =-J J drdY'[^ K (r,r')] 2 W(r,r'). 

If we define the functions /(r) = ^{^^(r)] 2 + 
z 2iS [Xb(j)] 2 } and g(r) = e lS xa(j)xb{t), then A can be 
written as 



A = A 



AA 



*AB , 



with 



A AA = - J d 3 rd 3 r'[f(r)rW(r,r')f(r'), 
Aab = I d 3 rd 3 r'[g{v)]*W{v,v')g{v') 



Aaa is strictly positive and Aab is strictly negative. 
Let the integrated charge of the product of like basis 
functions be defined as Qaa = J d 3 r{xA(r)] 2 r and for 
different basis functions Qab = / d 3 rxA( r )XB( r )- Since 
the product of two gaussians with exponent a centered 
on two neighboring atoms a distance d apart is a gaussian 
centered on the midpoint of the two atoms with factor 
e -ad li anc j there are three nearest neighbors, then we 
have, approximately, Qab ~ 3e~ ad ~/ 2 QAA- 

Since the dominant exponent for our calculations is 0.2 
(a.u.)" 2 , the ratio is ie~ aS I 2 = 1.4 > 1.0. The value of 
0.2 (a.u.) -2 for the exponent a does not depend on the 
chirality and size of the tube since it is a characteristic 
only of the local, spatial character of the bonds. There- 
fore, it should be generally true that the term Aab is 
dominant over Aaa and, consequently, A is negative. 
We expect the dark state to be generally lower than the 
bright state for carbon nanotubes, as we see in our ab 
initio calculations for several cases. This is so, because 
the state v + , which has the higher eigenvalue, E+, is an 
antisymmetric combination of the transitions at K and 
K' (the relative sign is given by the factor A*/|A|). 

For sufficiently long radiative lifetimes, exciton relax- 
ation to triplet states may also become relevant 01 ■ We 
therefore briefly discuss the splitting between singlet and 
dipolc-forbiddcn spin-triplet states. 

To calculate such splitting, we perform two BS calcu- 
lations with two different hamiltonians, one for the sin- 
glet state, H s = D + 2X — W, and one for the triplet 



(n, m) 


d (A) 


E\ u ~ &lg 


TTlS TJiT 

&lg ~ Ulg 


&lu ~ E lu 


(4,2) 


4.14 


0.090 


0.035 


0.345 


(6,4) 


6.83 


0.040 


0.010 


0.130 


(8,4) 


8.29 


0.025 


0.005 


0.075 



TABLE I: Energy splittings in eV. The third, fourth, and fifth 
columns are the splittings for the u — g singlet states, singlet- 
triplet for the gerade states and the single-triplet splitting for 
the ungerade states, respectively. 



state, given by H t — D — W. These hamiltonians differ 
only in the exchange term. Therefore, unlike the dark- 
bright splitting of the singlet state, which is due to the 
direct Coulomb interaction W, the singlet-triplet split- 
ting is due to the indirect exchange Coulomb interaction. 
Our ab initio calculations show that for the triplet states, 
the even exciton lies above the odd one [i.e., we have 
Ej > Ef u ], which is a trend opposite to that found for 
the singlet states. This behavior has been found also 
for zig-zag tubes |l6| . The above demonstration for sin- 
glet states does not apply to the triplet states. This 
is so because we have found numerically that the sim- 
ple four-band picture (two pairs of conduction-valence 
bands) does not hold. In fact, while for the singlet state 
99% of the weight is concentrated on the K and K', for 
the triplet states at least 6% of the weight is concentrated 
elsewhere. By evaluating A due to these other contribu- 
tions, we find a positive contribution which is enough to 
reverse the sign of A, forcing the even triplet exciton to 
lie above the odd one. 

The above symmetry analysis of the singlet states sug- 
gests new possibilities for increasing the quantum effi- 
ciency through symmetry breaking. It can be deduced 
from the above discussion and perturbation theory that 
the dark state becomes optically allowed, if the interac- 
tion which breaks the symmetry is of the order of magni- 
tude of the splitting E u ~E g m 0.1 cV. One way to break 
the symmetry involves applying external fields. How- 
ever, their required magnitude may be too large for re- 
alistic technological applications, as shown in Ref. [l7j . 
where magnetic fields are used. An alternative strategy 
could involve the use of interacting semiconductor nan- 
otubes. This is suggested by recent theoretical work on 
excitons in polimcric crystals, where interactions between 
different polimeric chains were shown to control the split- 
ting between dark and bright excitons mill. Wc have 
therefore reason to believe that bundles or arrays of semi- 
conducting nanotubes could be examples of systems with 
increased quantum efficiency and possible candidates of 
interest for practical devices. 

In conclusion, we have provided theoretical and nu- 
merical proof, using an ab initio, many-body approach, 
and separately, a four-band model, that dark cxcitonic 
states in semiconducting single- walled nanotubes exist, 



4 



that they have even Ru symmetry, and that they are 
always below the corresponding bright state. Further- 
more, the dark-bright energy for tubes in the diameter 
range of 0.4 - 0.8 nm is some tens of meV with decreasing 
value as the tube radius increases. We show that these 
results come from a four-band picture and that the split- 
ting is due to the direct Coulomb interaction between the 
electron-hole and the matrix element of the inter-valley 
interaction. We also calculate the singlet-triplet splitting. 
Here we find results of the same order of magnitude and 
the same decreasing trend with increasing radius as in the 
case of the dark-bright splitting. Finally, we have pro- 
posed, on the basis of the symmetry characteristics of the 
dark state, coupled semiconductor nanotubes as systems 
where the symmetry breaking induced by the intertube 
interaction leads to increased quantum efficiency. 
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